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Abstract 


A Power Bipolar Junction Transistor (PBJT) model is developed giving both 
static and dynamic behavior. The ordinary differential equations are derived based 
on the semiconductor physics. These equations take into account the recombina- 
tion, charging and discharging process in the lightly doped collector drift region. 
The model is implemented along with associated circuitry, like base drive, load 
(resistive and inductive), by writing a C— program. The static behavior giving 
quasi-saturation characteristics can be simulated with this model. The transient 
behavior of the model is studied with different base drives, time steps and error 
margins. Also charge profile in the drift region is presented at different instants 
during turn on/off. Experimental results , wherever available, are plotted along 
with simulated results. The matching between the two validates the model devel- 
oped. Also qualitative agreement with published work exists for various simulated 
results presented in this thesis work. 
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Chapter 1 


INTRODUCTION 


The Power Bii)olar Junction Transistor (PBJT) differs in many respects in contrast 
to tlu^ low-volt:i;;e BJ'I'. In PBJT there is a need for large blocking voltage in the 
off-stal.e and a liigh current carrying capability in the on-state. This leads to a 
dilfercnce in sta ucture and doping profile of a PBJT. Due to different structure and 
doping profile iJie I V characteristics and switching behavior of a PBJT is quite 
different from a low-voltage BJT. Due to above mentioned reasons low-voltage 
BJT models arc not suitable or applicable for the PBJTs. In this work a PBJT 
model is devcloi)ed which gives both static and dynamic behavior under different 
bias conditions. This model is derived and implemented from the study of various 
works publislu'd in this field as given in references [l]-[4]. 


1.1 Structure of PBJT 

There are many structures of PBJT e.g. n'^p n"^ , n^p u n"*" , n"*p v nn^ etc. The 
PBJT structure considered for modelling in this work is n^p v structure. The 
structure of the PBJT is shown in Fig. 1.1. 

The PBJT has a vertically oriented four layer structure . The vertical structure 
is preferred for power transistors because it maximizes the cross-sectional area 
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Figure 1.1: Structure of PBJT 


through which the current in the device is flowing. This minimizes the on-state 
resistance and thus tlie power dissipation in the transistor. 

Tlie doping levels in each of the layers and the thickness of the layers have a 
siguiflcant efl'ect on the characteristics of the device. The doping in the emitter 
layer is quite large (typically 10^® cm“®), whereas base doping is moderate (10^® 
ciu“'^). The /^-region called the drift region has a light doping level (10^'* cm“^). 
The 71''' region that terminates the drift region has a doping level similar to that 
found in the emitter. The thickness of the drift region determines the breakdown 
voltage of the transi.stor and thus can range from tens to hundreds of microns in 
extent. 

Practical ])owor transistors have their emitter and bases interleaved as narrow 
fingers to reduce the effect of current crowding, a phenomenon that can lead to 
device failure. This type of structure also reduces the parasitic ohmic resistance in 
the base current path, which helps to reduce power dissipation in the transistor. 
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1.2 Basic Model of PBJT 


Tiio PBJT cliaractc'ristics are dominated by the charge in the lightly doped collector 
layer (drift region) which is negligible or non-existent in a low-voltage BJT. This 
nmst be eoiisidcned and described in a generic PBJT model. The basic model is as 
shown in the Fig. 1.2 [5]. 



Mgure 1.2: PBJT Model 


The model consists of a standard BJT model and two additional elements, the 
collector resistance Rcc and a current source Icp- Rcc is added to take care of 
the voltage drop in the drift region and Icp basically describes the recombination 
taking place in collector layer. Icp is hole current from base into the collector. The 
recombination in rollcctor layer is not negligible in spite of high lifetime and thus 
must be incorporated in the model. These two elements along with standard BJT 
model are able (,o describe both static and dynamic characteristics of PBJT. 


3 




1.3 Basic Approach for Computer 
Simulation 

A widely ;\cc('pted a])proach to the fundamental modelling of semiconductor de- 
vices involves t he solut ion of the semiconductor device equations, which comprise 
Poisson’s efination and the hole and electron continuity equation. This, however, is 
often hnpraetieal as it takes a large CPU time and memory. In this work, we treat 
the PUJT wit.h -a, circuit model instead, and solve the circuit equations together 
with the nioch'l cciuations. The Newton- Ilaphson iteration method is used to solve 
the set of nonlinear differential equations. A second order integration method is 
used for implicit time stepping i.e. Gear-Shichman integration method. The con- 
vergence to a solution is controlled by keeping a very low error margin of about 
10~®. 'Phe basic algorithm flow chart is given in Fig. 1.3. 

First the circuit to be simulated, i.e., PBJT and the external elements like supply 
voltages, load, etc. are converted into a set of equations in the form of matrices. 
The DC analysis of the circuit is carried out to find the bias point. The matrices 
are solved using gaussian elimination method. For transient analysis the results of 
DC analysis sui)pty an initial guess. All the derivatives are removed using Gear- 
Shichman int('gration method. Time step is set to be about one hundredth of 
minority carrier lifetime. The equations are then solved at each time step using 
Nowf.on-Paphson iterations and error is worked out. Convergence is only achieved 
if error goes below 10"® at each node. The results thus achieved supply guess for 
next time step. 

The results thus achieved are checked and compared with the published results. 
Quite good agHHuneut is achieved in this. 
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Figure 1.3: Flow Chart for PBJT Simulation (Static and Dynamic) 





























Chapter 2 


MODEL DESCRIPTION 


2.1 The Static Model 


The basic equations governing the transport of carriers are the Poisson, the current 
density, and the continuity equations. As the collector layer is at least partially 
high injection in the on state {n^p), the following equations describe the minority 
carri('r distiibul.ion p{x) in the high injection region for a one dimensional current 
transport, which is assumed for the collector layer [5]. 


3p — 


b+1^ 



( 2 . 1 ) 


^ ^ V Op ^ P 

e dx dt T dt 


( 2 . 2 ) 


where e is the electron charge, b the ratio between electron and hole mobility, and 
D the ambipolar diffusion constant, j and jp are the total current and hole current 
density, respectively, 9^ is the recombination rate. A constant lifetime r is assumed. 

At steady state {d/dt = 0), eqns. (2.1) and (2.2) are reduced to an ordinary 
differential eejuation, given by 


= JL- = JL 

dx^ Dt 


(2.3) 
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where L is the ainbipolar diffusion length. This is a second order differential equa- 
tion. The .solution of this equation is a hyperbolic; function. 

p{x) = C 1 cosh j + C 2 sinh y (2.4) 

I-j Lj 

where C\ and C 2 are constants which can be found using boundary conditions 
[5]. Figure 2.1 gives the charge profile in the drift region of PBJT. The boundary 



Figure 2.1; Charge Profile in Drift Region of PBJT (a) Quasi-Saturation (b) Hard- 
Saturation 


conditions are given by the following equations. 


p{xi) = 7ii exp 



(2.5) 


dp _ ie + (6+ l)/cp 

r.{b + l)DA 


( 2 . 6 ) 


dp 

dx 






e{b + l)DA 


(2.7) 


where p{xi) is the carrier concentration at the base-collector boundary, is the 
intrinsic carrier density, Vt is the thermal voltage, and Vbc is the internal base- 
collector voltage, i.e., between the base and collector of the standard transistor of 
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Fig. 1 . 2 . Xi and X 3 are the edges of the injected charge towards the base and the 
n"*" colh'ctor contact, respectively as depicted in Fig. 2 . 1 . is the collector current, 
Icp is the recombination current and A is the collector base junction area. 

If the charge distribution is such that after turn-on ^3 is not able cover the 
entire width of the collector drift region then the transistor is operating in the 
quasi-saturation region; else in the hard-saturation region as depicted in Fig. 2.1, 

Thus there are two possible solutions for the collector layer. Now if W is the 
width of the collector drift region then the solution for quasi-saturation are as given 
below. 

Fl'om eqns. (2.4) and (2.5) 


C'i=p(:ri) 


( 2 . 8 ) 


where x = xi =0 

Also for Quasi-saturation pix^) = 0 . Thus from eqn. (2.4) 
p{xi) cosh ^ 


C2 = 


sinh ^ 


Substituting Ci and C 2 in eqn. (2.4) , we get 
sinh 23^ 

Applying boundary condition of eqn. (2.7) on eqn. (2.10), we get 


Xz = L sinh 


-1 


e (6 -f 1) DAp{xi) 

%c T 


< W 


(2,9) 


(2,10) 


( 2 , 11 ) 


Applying boundary condition given by eqn. ( 2 . 6 ) on eqn. ( 2 . 10 ), we get 
ePApjxi) cosh (f) -1 


sinh ^ 


( 2 . 12 ) 


Rcc consists of two parts, i.e., resistance of the charge induced part(i?ci) and re- 
sistance of the unmodulated part(i?c 2 ) of drift region. Now Rd and Rc 2 can be cal- 
culated by cahndating the voltage drop in the respective drift region parts(V'ci, Vc 2 ). 
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If the (Iroj) (hi(' to the recombination current Icp is neglected then Rd, i?c2 can be 

cmIcuImI t'<l in (In' lolluwiuf.’ iiumiu’i’ [f)], 




111 

■>c 


Irln 



Nd J 


Rc2 — 



Rro 



TI-" 


Where Rco is tlio resistance of drift region without injection of any charge and No 
is the doping density of the drift region. Rcc is given by the following relation. 


Rrc — Rcl + Rr2 


substituting values of Rd, Rc 2 , we get the following relation for Rcc- 



and the solutions for hard-saturation are 


(2.13) 


X 3 — w 


(2.14) 


Cl and C 2 can again be evaluated using eqns. (2.5) and (2.4). Ci is same eis in 
eqn. (2.8) and C-z is as follows. 


Co 


p(a' 3 ) — p{xi) cosh 

Ij 


w 


(2.15) 


p{x 3 ) can be evaluated using boundary condition given by eqn. (2.7) and the 
value thus obtained can be used to evaluate Cz- Substituting Ci and C 2 in eqn. 
(2.4) , we get 


p{x) = 


p{xi) cosh 


W-x 


— ISLld. SlTlll ~ 

L e(6+l)PA I. 


cosh ^ 

iJ 


(2.16) 


Applying boundary condition given by eqn. (2.6) on eqn. (2.16) , we get 


, eC.4p(.T,) ^ ic cosh(f)-l 

L L 6 + 1 coshf 


(2.17) 
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(2.18) 


]±, f pCrO+Tfp ] 

VP(-'^3) + A^i)/ 


Again jR-rc is calculatod in the similar way as in the quasi-saturation case, only 
difference being that, now entire drift region is charge injected so there is only one 
type of voltage drop which can be calculated knowing p(xi), p(x 3 ) as given in [9]. 
In calculation of I?cc voltage drop due to recombination is neglected as it is very 
small as compared to other voltage drops in the PBJT. 

The circuit model for modelling the static characteristics of PBJT is given in 
Fig. 2.2. The ecjuivalent circuit of PBJT model of Fig. 2.2 is given in Fig. 2.3. 
The standard transistor of basic model is represented by the Ebers-Molls equivalent 
circuit. 


node c 



Figure 2.2: PBJT Circuit (Static Case) 


A constant current source It is assumed. From the figure it is clear that there 
are three nodes (5, z, c). The node b is the base node, node i is the internal collector 
node, i.e., collector node of standard transistor, and node c is the collector node ' 
of PBJT. 1 ce is variable voltage source. For nodes b and i, KCL is written and 
Branch Constitutive Equation (BCE) is written for node c. 


10 





Figure 2.3: PBJT Equivalent Circuit (Static Case) 

Mvb, Vi, ic) = he (1 - cxf) + he (1 -ar)+lcp-h = 0 (2.19) 

^ 2(^61 ‘^ii *c) = Icp he he ~ 0 ( 2 . 20 ) 

— l^ce -^cc ~ 9 (2.21) 

where Icp , /?,c iue given by equations (2.12) or (2.17), and (2.13) or (2.18) respec- 
tively as the case may be and he and he arc as given below. 

he = h [oxp - l] ( 2 . 22 ) 

he = h [exp - l] (2.23) 

where h is the reverse saturation current and /i,/ 2,/3 are the functions of three 
unkjiowns. 
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The unknowns for the circuit simulation are node voltages at node b and i and 
node current through nodes c and i i.e. ic, the collector current. An initial guess 
is assumed for the unknowns and equations are solved using Newton Raphson 
iteration techniciue as given below. 

The function values and the partial derivatives of each function with respect to 
the unknowns are evaluated to form the following equations in matrices form. 


^ SA iA 

dl% dl'r 


dj2 ^ 

dvi dvi die 


dfs Oh Oh 

Ovi, OVr Otc 



1 

> 

1 



Aui 

= — 


Aic 



Mvlvlil) 


(2.24) 


where Uj, vf, are the assumed values of Vb, Vi, ic , first matrix is called the 
Jacobian and the column vector on left is the error vector. Equation (2.24) is 
solved using Gaussian elimination technique and value of the errors are worked 
out. If error for all the three unknowns is within desired error margin then the 
process stops else a new guess is generated by adding the error for each unknown 
into the previously assumed value. The iteration process continues till error matrix 
is reduced to within limits of desired error margin. This technique of solving non 
linear equations is known as Newton-Raphson Method. Vee is varied to get the IV 
curve for a particular value of base current h ■ This process is applied for each 
new value of 1 cc and the initial guess for each case is the solution at previous Vee 
value. 


2.2 The Dynamic Model 

The circuit in Fig. 2.2 describes also the dynamic behavior, but some variation 
must be considered. The formulation of Rcc remains approximately valid. Icp now 


12 



contains the variation of the charge in the collector layer. 


■frp — 


Q dQ 
T (it 


(2.25) 


whore Q is t.lu' ston'd charge in the lightly doped collector layer given by 


O ^ c /I 



(2.26) 


To describe p(x), eqns. (2.1) and (2.2) are reduced to a single equation. 


d'^p If dp\ 
dx^ ~ Dt dtj 


(2.27) 


The above equation must be reduced to an ordinary differential equation to get 
a network model. The approximation made to realize this model is [5]: 


dp _ p 

Tt~T 


(2.28) 


where T is a function of time but not of position. Thus eqn. (2.27) can be written as 
given below. This approximation is made for the simplicity of the model. Comment 
will be made on this approximation in chapter 4. 


d'^p 

dx"^ 




(2.29) 


where A is the newly introduced effective ambipolar diffusion length. All the bound- 
ary conditions given by eqns. (2.5), (2.6) and (2.7) remain valid. Equation (2.29) 
can now be solved analytically. 

In the st.atic cfise, 


It ~ dt 


0 ^ = 0 


(2.30) 


A is identical with L. The solution forp(a:) is same as given by equations (2.10) or 
(2.16) as the case may be. 
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For turning on of PBJT: 


(IQ 

(it 


> 0=^^>0 

dt 


> 0 


(2.31) 


A is smaller than L, now p{x) will be a hyperbolic function as in the static case 
where A will replace L, i.e, eqns (2.10) to (2.18) will have L replaced by A. 

For turning off of PBJT {dQ/dt < 0) there are three different cases: 


(a) Charges are still injected into the collector layer {Qp > 0), 


dQ ^ Q dp P ^ 1 

dt T dt T T 


(2.32) 


A is now greater than L and p{x) is still a hyperbolic function where A replaces 
L as for turn on. 


(b) No charge is injected into or extracted from collector layer {Qp = 0). 


dQ 

dt 


Q dp 

— -f 

r dt 




(2.33) 


p{x) is now a straight line given by 


p{x) = p{xi) - 


{x - Xi) 


e{b+l)DA 

(c) Charge is removed from collector layer (Jcp < 0). 


(2.34) 


dQ Q dp p 
^ < “7 


r 


< -1 


(2.35) 


Now the sign of eqn. (2.29) is negative and p{x) is a sinusoidal function given 
by 


p(:;:,t) = CiCOS Q + ^2) 


(2.36) 


again Ci , C 2 can be determined by applying boundary conditions mentioned 
before. Also Icp can be evaluated in similar manner. Finally solution is: 


cos 

sm [f) 


re - (x3- I a) 
A 


(2.37) 
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(2.38) 


‘-cp 


6 + 1 


1 — cos 


fa' 

A. 


where 


.r:i — A si if 


p{x\)e (6 +1) D A 


XL 


(2.39) 


E(niation (2.28) is an approximation and is assumed only for a certain region. 
The model would be exact if the whole collector layer is divided in many such 
regions, but this api)roach will be very CPU time consuming. Therefore number of 
divided regions must be low to reduce the complexity of the model, i.e., the number 
of regions (each with its own time constant) should be as few as possible. Usually 
one time constant for the whole high injected region is sufficient. During turn off 
the reverse base current is sometimes very high so that [L + (6 + l)Icp] < 0. This 
means dp/dx at a: = xj is positive , the sign of slope of p{x) changes in the high 
injection rc'gion, and p{x) has a cusp [3]. This drive is called a strong drive. Here 
only one T (sinusoidal function) for p(x) would violate the balance of electron and 
hole current and characterize the process of discharging inaccurately. In this case, 
two sinusoidal functions (regions) are necessary to describe the charge distribution 
from the cusj) towards the base and the n'*' collector contact. 


P2 cos for X\ <X <X2 

P 2 cos for X2 < X < X3 , 

where p 2 is i lu' minority carrier concentration at the cusp and X 2 is the distance of 
the cusp from base collector junction. Also 1^, is evaluated in similar manner as 
described previousl}^ 




_ + A 2 

~ e{b+l)DA 


(2.41) 


1 


cp 


6+1 




(2.42) 


The value of A can be worked out by the solution of eqn. (2.25) which is given . 
by 


Q — T Icp (^1 — + Qo 6 ‘ 
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(2.43) 



where Qo is the charge in the drift region at t = 0. Now the value of Q can be 
worl<('il onl iisiiif.'. (’i|n {? 2G), lly Hiili«til.iitiiip, Uiiri v-'iliic' (>[ Q in cqii, (2.4G) arirl 

solving for A, we get 

A = L''-’ (l - (2.44) 

For t.iirning on Qo = 0 and for turning off Qo = t Qpo- Here Icpo is the value of 
Icp obtained after turn on is complete. Also if there are two time constants, we get 
the following relation. 


Ai A2 = L" (l - + 


QoDe-^l^ 

Icp 


Also by the condition at the cusp [3]. 


In = hl, 


In — V Icp 


AQi 

At 



AQ2 

At 


(2.45) 

(2.46) 

(2.47) 

(2.48) 


where /„ is the electron current, /p is the hole current, Qi and Q 2 are the charges 
on both the sides of cusp. The above equation set yields the value of Ai, Ai. Thus 
all variables for transient analysis are known. 

The circuit model for transient simulations for PBJT takes into account a base 
drive supply voltage and also a load at the output node, which can be resistive 
or inductive as the case may be. Also a snubber capacitor is added between the 
collector and ('initter of PBJT to fasten the switch off time. In practical circuit 
usually a diode is used instead of a capacitor. The circuit model is illustrated in 
Fig. 2.4. The base drive voltage is positive during turn on and polarity changes 
instantly during the turn off. 

To drive the ecpiivalent circuit of Fig. 2.4, each diode of Fig. 2.3 is modelled as 
shown in Fig. 2.5 [G]. 
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Figure 2.5: Diode Equivalent Circuit 

where 

/,i = /<i/ + /jc + /r (2.49) 

i.e., the total current flowing through diode constitutes of three subcomponents 
which can be modelled as 
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(a) DC current component I^f given by 


l,,f = Js (exp ^ - i) 


(2.50) 


(b) Junction capacitance current component Ijr, given by 

r 


(2.51) 


dt 

where Qjc is the charge associated with the junction capacitance Cjc given 
by 

C, 


C,c = 


']0 


( 1 -^) 

Qjc ~ ^ CjgViji ( 1 


VhJ 


(2.52) 


(2.53) 


Here Cj„ and vi,i are zero bias junction capacitance and built in junction 
potential resi)ectively. Also Vj is the applied potential. 


(c) Diffusion current component Ir, given by 


Jr = T 


didf 

dt 


(2.54) 


Thus both the diodes of Fig. 2.3 can be expanded to generate equivalent circuit 
of Fig. 2.4. The equivalent circuit is given in Fig. 2.6. 

The procedure for carrying out transient analysis again consists of writing the 
nodal equations and then solving them. Clearly there are three nodes as indicated 
in Fig. 2.G. For each node KCL and BCE are written. As we are concerned with 
the the transient analysis the equations will contain (d/dt) terms as well. These 
terms are apijroximated using Gear-Shichman approximation technique and then 
equations are solved in a similar manner as in the DC analysis case. However 
the choice of time-step is very important because a bigger time step can result in 
loss of transient behavior and a very small time step can result in too much of 
computation time. Usually time step is fixed at one hundredth of the minority 
carrier lifetime. The following set of equations illustrates the full procedure. Two 
cases are being considered, one for resistive load and other for inductive load. 
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Figure 2.6: PBJT Equivalent Circuit (Dynamic Case) 

2.2.1 Resistive Load 

Writing KCL and BCE for each node. 

Node b : 

f\ — ^be Ijc “h Ide ^r^bc ^b ^bc d” ^cp d" Ijc d" ide ^bc “ ^ (2.55) 

where hr., he fhe dc current components, Ijc, Ijc are the junction capacitance 
current components, and he, Ide are the diffusion capacitance current components. 

^ /l = Ibc.n+i d- Qje„+1 d" rie he^+i ~~ ^bcn+i ~ Ibn.+ i d" hcn+i d" ^cpn+l 

d"Qjc„+i ‘^bchc„+i ~ Iben+i ~ ^ (2.56) 


The Gear-Shichman approximation is given by 


V' 


^ Y —X 


(2.57) 
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where ri+ 1 denotes the value of X at n + 1 time instant, h is the time step, and X' 
(lenott's iliM'ivativc of A’ witli reHj)(‘ct to time. All the derivatives of crpi. (2. GO) can 
be approximated using eqn. (2.57) technique and formulae given by eqns. (2.52) 
to (2.54). The final nodal equation thus obtained is given ahead. 


/i — hc„+i 


3 Thf, 

Til 


+ 1 - a/) + + 1 - or.) - ^ 

‘^■nc 


4 -— T 
+ 2 /, Wi 






L Qjsn 2 h ^ o h 


n+l d" ^Cpn+1 


A 

2h 


_2 _1_ 


(2.58) 


Similarly at Node i 


/a -O'/ ( 2/i +l) + |^ 2/1^'"'"-* 2/1^^'='*+* 


'f' j^QjCn ^cpn+l *C„ + 1 —0 (2.59) 


Node c 


h 


i 


C'n + 1 


r 3(7 

2 h 



C 

2/i^' 


= 0 


(2.60) 


BCE: 


^4 — Vf,5 ^6n+l -^frn+l ^ 


(2.61) 


h = ’^Cn + l — ^»n+I ~ *Cn+l -^CC ~ ^ ( 2 . 62 ) 

/g = brr — t.'m+i ~ — 0 (2.63) 

h = Uc„+i " “ 3 ^‘=’‘ 3 ~ ^ 

The above equations (/i to fj) represent the complete set of equations required 
to achieve the transient solution of the PBJT circuit. The Jacobian can be worked 
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out, aiui function values can be calculated around the initial guess provided by the 
DC analysis. For the DC analysis the capacitor will be open circuited and load 
resistance value can be added to Rcc- The procedure for solving the equations re- 
mains same, i.e., involving Newton-Raphson iteration method as explained earlier. 
At each time step solution is thus obtained for transient response. 


2.2.2 Inductive Load 

The equation set can similarly be generated for the inductive load. In this case 
load is not purely inductive, but a resistive load is connected in series with inductor 
for proper biasing of the transistor as shown in Fig. 2.7. 



O 6 


Figure 2.7: Inductive Load Circuit 


The changes will be in the function values of /a, fe and due to the addition of 
node vi two more functions will be evaluated corresponding to unknowns vi,Iri. 
These function values can be obtained with the help of the following equation. 

h.„ = (Kc - + pL,-l h.-. (2.65) 

This value of will replace /jr, in the equations of to give the set of 
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equations for transient response with inductive load. During DC analysis inductor 
will be short circuited. The final set of equations is as given below. 




( 2 . 66 ) 

( 2 . 67 ) 

( 2 . 68 ) 

( 2 . 69 ) 


The remaining function equations i.e., fi, f2, fi, 15 , fr remain same. Remaining 
procedure is identical with the resistive load case. In this way the transient analysis 
can be performed for PBJT for various types of loads. 
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Chapter 3 


RESULTS AND COMPARISON 


The model of PBJT described in the last chapter was implemented with a C- 
program. The various model parameters include the standard low power transistor 
model parameters and five additional parameters viz. W, A, r, N^, and Rco- All the 
above parameters are already described in detail in the last chapter. The values of 
the various parameters were obtained from [1],[2]. Using all these values the results 
for static and dynamic behavior were derived by running model program. 

The model parameter values are listed in Table 3.1. For parameters not given 
in [1],[2] some reasonable values were chosen. 


Table 3.1: Device Parameters for Static and Dynamic Analysis 


DEVICE PARAMETERS 

Temp = 24° C 

hpBO 

Of 

Oer 

Nd 

W 

D 

m 

r 

Mn 

I- 

Cj„ 


'Hjc 

Rco 

B 

Vbi 




10*^ 

cm~® 

fj.m 

cm^/s 

D 

MS 

crn^/Vscc 

pA 

nF 

MS 

MS 

n 

1 

V 

16,5 

0.94 

0,5 

1.16 

50 

22.3 

0.174 

15 

1300 

2.0 

0.46 

0.01 

3.0 

2.69 

3 

0.8 


Note that 17 is the ambipolar diffusion constant and r is the miniority carrier 
lifetime in the lightly doped collector layer. 
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3.1 DC Characteristics 


The DC characteristics were obtained for the model PBJT as mentioned earlier 
(given in Fig. 2.2). These curves along with the results from [1],[2] are given in 
Fig. 3.1 and Fig. 3.2. The curves are obtained for different base drives as indicated 
in the graphs. 



Figure 3.1: PBJT DC Characteristics for Higher Value of 4 (a) Model Results (b) 
Results of [1] (c) Results of [2] 


The curves of [1] do not take into account the collector lifetime r, whereas effects 
of r have been accounted for in [2]. The model results are in good agreement with 
results of [2], which further shows the validity of the model. There are two distinct 
regions in the eharacteristics of Fig 3.1 and throe regions in Fig. 3.2. The first 
slope indicates hard-saturation region, the second slope indicates quasi-saturation 
region, and the final one being active region. 


In the quasi-saturation regime, the behavior of the PBJT is basically dominated ' 
by the charge stored in the collector layer. Also the collector layer lifetime is very 
important in determining the PBJT behavior because recombination current Icp 
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Figure 3.2: PBJT DC Characteristics for Lower Value of h (a) Model Results (b) 
Results of [Ij (c) Results of [2] 

is affected by this lifetime. All these factors are taken into consideration in the 
derivation of this model. 


3.2 Transient Behavior 

3.2.1 Resistive Load 


The transient characteristics are also obtained by using the procedure mentione^j 
in the last chapter. The various model parameters and the component values are 
taken from [l],[2] for comparison of results. The PBJT is switched on and off by 
J'-Pplying a base drive of +2 V and -2 V respectively. The change in the base driy^ 
voltage is assumed to be instantaneous, i.e., from 0 to +2 V for turning on an^i 
from +2 F to ~2 V for turning off. The PBJT is off at time t = 0“ and turn 
is initiated at time t = 0. Figures. 3.3, 3.4, and 3.5 give the various current an^} ' 
voltage waveforms obtained by simulation of the circuit shown in Fig. 2.4. Her^ 




A n Cf A 



base drive' change's polarity at t = 20Atsec during turn off. The arrow in the time 
axis iiidie-ate's (he lime at which turn olf is initiated. 



Figure 3.d: Transient Current/ Voltage Waveforms for Ri = 3.6 Q, {Ic, Fee) 



Figure 3.1: Transient Current/Voltage Waveforms for Ri = 3.6 (/&, Vbe) 


All the wave'le)rms are obtained for the purely resistive load of 3.6 Vec of 3.7 V 
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FiRiirc 3.5: Transient Voltage Waveforms for Ri = 3.6 Q, {Vce,Vci) 


and snublx'r capac it.or of 1 iiF. The base biasing resistor is of 5 ohms. A saturation 
resistance of 0.085 ohms is assumed for plotting Fee waveform. Figure 3.5 shows 
the comparison l)(‘two(>n the internal and external collector voltage. The time step 
of r/100 is talum for simulation of currcnt/voltagc transients. The error margin 
was fixed at lO"-’ for calculation of all the currents and voltages. The effects of 
changing the error margin, base drive, and time step are dealt in detail later in the 


chapter. 

Wo now con.paro llu- waveforms of Fig. 3.3 with the experimental results taken 
from [1]. Figure 3.6 gives the results tor the turning on of PBJT. In [1], the turning 
on process was modelled assuming two distinet time constants. The model results 
are in quite good agreement with results of (I). Both the curves are obtained for 
load rcsistimce value of 3.6 SI. The compatisou curves validate the model results 


for the turning on process. 

In 11] the turn on transients were given for another value of load resistance 
as well. The simulated results are aUo worked out for this case where the load 
resistauce is of 20.8 S! and K. of 21 V. In Fig. 3.7, the comparison between the 
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Figure 3.7: Comparison Curves with Ei = 20.8 Q (a) Simulated (b) Results from 
[ 1 ] 


Hero also there is good agreement between the simulated and experimental re- 
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suits. This fm tTor validates the model. All waveforms for a resistive load of 20.8 
ohms are given in h'igs. 3.8 to 3.10. 



Figure 3.8: 'Fransient Current/Voltage Waveforms for lii = 20.8 (/c, Kc) 



Figure 3.9: Transient Current/Voltage Waveforms for Ri = 20.8 Q (4, He) 
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Figiue 3.10: IVansient Voltage Waveforms for Ri = 20.8 (V^, Ki) 


The turn olf result s cannot be validated due to the non-availability of the exper- 
imental data for the model transistor. However, a good (jiialitativc agreement is 
achieved by comparison of the shape of the various waveforms with the usual shape 
of turning of rt'sults of PBJT as given in [7]. The model was tested under various 
bias conditions with different combinations of supply voltage and loads. During 
the various simulations it was noted that a proper selection of bias point is crucial. 

3.2.2 Inductive Load 

The model transistor’s transient characteristics were also checked with inductive 
load (Fig. 2.4 and 2.7). All the model parameters remain same except the snubber 
capacitor in this <‘asc was taken to be 10 tiF. A resistive load of 24 D was connected 
in series with the inductive load for proper biasing of transistor. The supply voltage 
Vcc wiis kept at 25 The same base drive was used to switch the transistor on 
and off. Again switching off was initiated at time t = 25 /j,sec. The time step and 
the error margin setting remains same as in resistive load case. Figures 3.11 to 3.14 
give the various current and voltage waveforms for an inductive load of 10 ixH . 
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Figure 3.11: 'IVan-sient Curront/Voltage Waveforms for Inductive Load (/c, K*) 



Figure 3.12: Transient Curreiit/Voltage Waveforms for Inductive Load (/6,Hc) 


The period of the oscillations during the turning off of PBJT in Fig. 3.11 is a 
function of L and C. In practical cases, as the capacitor will be replaced by diode, 
oscillations will lujt occur. Also it is quite evident that turning on/olf takes longer 
time with the addition of inductive load as compared with purely resistive load. 
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Figure 3.13; Transient Current Waveforms for Inductive Load [IiJcc) 



Figun> 3.M: 'rransient Voltage Waveforms for Inductive Load (VccVd) 


Tlie waveforms of current through inductor and capacitor (//, Ice) are given in 
Fig. 3.13. During the turning on h follows h and /cc = 0 ) current flows 

through the snubber capacitor. During turning off when 4 goes to zero, all the 
current goes thiough the capacitor. This is very clear from Fig. 3.13. Basically 
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the KCL at. tlu^ node c is always maintained, i.e., /l = /c + Icc- 


'riu' eouipaiison Ix'tween tlie internal and external collector voltage (Kce, Vd) is 
given in d.M. d'ln' dilleK'iut' iu the voltag('.s is (lue to ilie voltage dropped in 
the drift region, i.e., across the resistance i?cc in the model. 


As l.here is no experinu'ntal data available in literature for transient behavior 
with inductive' load, we ('ould not validate the model quantitatively in this case. 
How(n-('r, t.he (pialitat ivo agreement exists by comparing the shape of waveforms 
with [7]. 


3.3 Effect of Varying Base Drive 

The transient lu'havior of the device depends significantly on the base drive. A high 
base current rc'ducos the switching losses for turn on but will increase the switching 
off time autl losses for the turning off. This is due to the fact that though higher 
base current will switch on the PBJT faster, during turn off the storage time as 
depicted by base' curnnit waveform will increase resulting in higher switching off 
time. H('.vcrsc base drive will also determine the switching off time. The PBJT can 
be turned off by making base supply voltage zero or negative. The turn off time is 
biisically function of the charge stored in the drift region because that much charge 
is to b(i r('mov(>d to turn off the transistor. The stronger the negative base drive 
less time will Ix' tiik(*n for transistor to turn off. 

The above effects were verified with the model developed. The model transistor 
was switched on by 1 going instantaneously from 0 to 2 V and in another case 
from 0 to 4 \ ' at time t=0. The reverse base drive was kept same for both the cases 
i.e. at -2 P. All the remaining parameters and component values were also kept 
constant. Figures 3.15 to 3.19 give the various comparison curves for the above 
mentionetl two c{iscs. 

Figure 3.15 gives the expanded view of the turning on process. Clearly by taking 
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Figure' 3.15: EiiVet of Different Base Drives on during turn on (Expanded View) 



Figure 3. 1C: Effect of Different Base Drives on Jc during turn on 

a strong base drive the process of turning on has become faster. The iV Vbb drive 
makes the transistor reach the 90 percent of Ic value faster than 2 V base drive. 
The effect of this decrease in turning on time on the turning off time can be seen 
in Figs. 3. 1C to 3.19. All of these waveforms correspond to the resistive load value 


34 






Figm t' .1.17: FliVct of Dilfcrcnt Biise Drives on h during turn on 



Tinic(//8ec) — )■ 

Figure ,1.18; IClfect of Different Bjisc Drives on V{,„ during turn on 

of 20.8 ohms. Also the time step is set at r/100 and error margin at 10“®. 

hVom nil th<' above waveforms it is clear that, a stronger forward base drive 
improvc's turn on !>ut dt'grades turn off performance. Thus the choice of base drive 
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Timc(/jsec) 


Figure 3.19: Effect of Different Base Drives on Vee during turn on 



Figure 3.20: Effect of Different Base Drives on Ic during turn off 


for turning on of PBJT is very critical. In the following analysis it is shown that 
how the turning off base drive affects the dynamic behavior of PBJT. 

The model transistor was switched off with three different base supply voltages, 
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Figure 3.21: Effect of Different B?ise Drives on h during turn off 



Figure 3.22: Effect of Different Base Drives on Kr during turn off 


i.e., 0 V, —2 V' —4 V. In all the three cases the transistor was switched on with 
same base drive, i.e., 2 V oi Vbb. All the remaining parameters were kept same. 
Figures 3.20 to 3.23 give the various current/voltage waveforms under different bias 
conditions. 
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Figure 3.23: Effect of Different Base Drives on Vbe during turn off 

It is evident from these figures that as we keep on increasing the reverse base 
drive the switching off time decreases. This is to be expected as a large reverse 
bias removes the excess charge faster. 


3.4 Effect of Changing Time-Step 

The choice of the time step is very crucial for transient analysis. If the time step 
is too large then the transient behavior of a device can be missed out and on the 
other hand if time step is too small then the number of time instants for which 
calculation of currents/ voltages is required will be too many, resulting in increase 
in simulation time. Thus there is a need to select optimum time step. 

The transient behavior of PBJT is governed mainly by the charge in the drift 
region. This charge is a function of collector drift region lifetime t. Thus if time- 
step is chosen less than the value of r the transient behavior can be modelled 
properly. On repeated simulations it was found that the value of time step of 
r/lOO is quite reasonable for PBJT transient analysis. The effect of varying the 
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time step on the transient behavior is given in Figures 3.24 to 3.27. 



Figure 3.24: Effect of Change in Time Step on Ic 



Figure 3.25: Effect of Change in Time Step on h 


The turn on transients were relatively unaffected by time step. Hence, only the 
turn off transients for various currents/ voltages are given in previous figures. It is 
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Figure 3.26; Effect of Change in Time Step on 



Figure 3.27: Effect of Change in Time Step on Vte 

interesting to see that beyond a certain time step, i.e., r/200 the solution does not 
change significantly. Thus time step smaller than this will only waste CPU time. 
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3.5 Effect of Changing Error Margin 


The selection of error margin is very important in circuit simulation. This is because 
the equation solving procedure is iterative in nature. The procedure is said to have 
converged if the variables differ from the previous iteration by less than the specified 
error margin. Thus convergence is a function of acceptable error margin. A high 
value of error margin will lead to faulty results and on the other hand a very low 
value will lead to a lot many iterations where sometimes oscillations can occur 
giving no result at all. Therefore, there is need to select optimum error margin 
such that number of iterations are less and the results are reasonably accurate. 
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Figure 3.28; Effect of changing Error Margin on Ic 


In the present model development the effect of changing error margin was studied 
to select the optimum value. Keeping the various model parameters, component 
values and supply voltages constant, the error margin was varied from 10~^ to 
10~®. This range of error margin was selected because various currents/ voltages in 
PBJT are much higher than this range, i.e., error possible with 1 amp/volt will be 
in mA/mV respectively. The waveforms of various currents/voltages thus obtained 
are given in Figs. 3.28 to 3.31. From the waveforms it is evident that beyond the 
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Figure 3.29: Effect of changing Error Margin on It, 



20 22 24 26 28 30 

Time(/Ltsec) — ¥ 


Figure 3.30: Effect of changing Error Margin on Vce 

error margin of 10"^ , the solution does not change noticeably. Thus the error 
margin value of 10“® is most suitable for the present circuit simulation. 

All the waveforms presented are for turn off transients only as the results in 
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Figure 3.31: Effect of changing Error Margin on V^e 

the turning on transients did not show significant changes for the selected range 
of error margin. Also it is interesting to note that there is a remarkable shift in 
solution when we change from error margin of 10“^ to 10~^ whereas further change 
in error margin results in very small shift in the solution and finally solution does 
not change appreciably beyond an error margin of 10~®. This was the reason for 
selecting this value for PBJT simulation under different bias conditions. 


3.6 Charge Profile During Turn on/ofF 

The dominant difference between a PBJT and a low voltage transistor arises due to 
the stored charge in the drift region. During the turn on process the charge keeps 
on building up till it reaches its steady state value as decided by the bias conditions 
and during turn off process the stored charge gets removed as per the reverse base 
drive condition. We have already seen how the time for removal/storage of charge 
varies with the t bias conditions. However, it is very interesting to actually capture 
the charge i)rofilc at different time instants during turning on/off. 
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Figure 3.32: Charge Profile in Drift Region of PBJT during turn on 



Figure 3.33: Charge Profile in Drift Region of PBJT during turn off 


The charge profile in the collector drift region was captured during turning on/off 
at the various time instants. These are shown in Figs. 3.32 and 3.33. All the curves 
correspond to the load of 20.8 and remaining parameter values same as before. 
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From both the figures it is very clear how the PBJT charge builds up during 
turning on and how it discharges during turning off. During removal of charge 
after the time instant 23.5 {xsec the depletion region forms at collector base junc- 
tion and charge profile changes to bisinusoidal function. This profile is there for 
very small time, therefore it was not captured.Thcsc two figures helps in actually 
understanding the buildup and removal of charge in the drift region of PBJT. 
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Chapter 4 


CONCLUSIONS AND FUTURE 
WORK 


The present model has been developed taking into account the recombination cur- 
rent and voltage drop in the lightly doped drift region. In the development of this 
model, several approximations were made which need to be discussed with respect 
to its application. 

The first approximation is that the device is one dimensional, and none of the 
two dimensional aspects was considered during model development. This was as- 
sumed to keep the governing differential equations simple. However, practically 
this assumption is valid only if the device is very fine structured, which is normally 
the tendency for the current development of the power devices. The devices with 
larger structure have two dimensional effects such as current crowding and these 
cannot be neglected. These are not incorporated in the present model. One pos- 
sible solution is to use a parallel connection of several PBJT models as shown in 
Fig. 4.1 [5]. 

Here the real transistor is separated into different slices. Now each slice can be 
described as a single PBJT model. All the collectors and emitters are connected as 
shown in the figure. Base resistors are inserted between the bases whose resistance 
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Figure 4.1: Two Dimensional PBJT Model 

values are controlled by the injection in base, the base-emitter and base-collector 
voltages of the two adjacent transistors. In this way 2-D effects in PBJT model 
can be implemented. 

The other important approximation is dp/dt = p/T for certain regions as dis- 
cussed in Chapter 2. There are four types of charge profiles possible during turn 
off. These are hyperbolic, linear, sinusoidal and bisinusoidal. During turning off 
the charge profile changes from one form to the other at different time instants. 
The change of its form depends on the bias and the carrier lifetime. The stronger 
the change, higher is the deviation of the approximation. Thus we can say that 
the inaccuracy of the model increases with the switching speed. However, in the 
absence of this approximation p is a function of position as well as time, and eqn. 
(2.27) cannot be reduced to an ordinary differential equation. Thus to keep sim- 
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j)licity ill tlie inodcl, approximation is made and as discussed in last chapter, results 
are qualitatively correct. 

Besides the stored charge in the collector layer, there are stored charges in 
the base and emitter as well. The space charges are described as the depletion 
capacitance. The excess charge in the emitter is negligible due to high doping 
whereas, the charge in the base has very complex characteristics due to structure, 
doping and the biases. This is due to the reason that base has no uniform doping 
and base current is two dimensional, one in the collector emitter direction and one 
in the cross direction. As the charge in the base is much lower than the charge in 
the collector drift region, the approximation in the standard model of the diffusion 
charge can be used. This was done while developing the equivalent circuit of the 
PBJT to give the transient behavior. 

We now briefly discuss parameter extraction which is very important for device 
simulation. The model contains a standard transistor model and two additional 
elements (Icp,Rcc)- Besides the parameter chosen, there are five additional param- 
eters for the two new elements {A,Nd,W,Rco,T). The first four parameters can 
be estimated from the device processing details. The most important parameter 
T can be found out using quasi-saturation analysis as given in [4]. The remaining 
parameters pertain to the standard transistor model and can be evaluated as given 
[5]. As the parameter extraction was not the main task in the present work, the 
detailed extraction procedure is not presented here. 

In the present work the PBJT is treated as a normal circuit element, i.e., its 
equations are solved along with the normal circuit elements like resistors, capacitors 
etc. [6]. The equation formulation for the PBJT is done by studying the physics of 
the device [7], [8], [9]. This particular procedure is not only valid for PBJT but can 
be extended for any low voltage or power device. Also this is the basic procedure 
adopted in any circuit simulator like SPICE. 

This model does not take into account the eff'ect of parameters changing with 
the operating temperature. Also the second breakdown effect is not included in 
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the model. These aspects can be exploited further in future work. As such this 
uiodt'l n good d(’Higuiug tool for power bipolar trauHistor having 

geometry. Also as a future work this model can be incorporated in SPICE as 
a separate model for PBJT. This can be done by adding the various equations 
prcseiitod for PBJT in the standard transistor model of SPICE. 
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